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Abstract 

The diffusion of a fractional Brownian particle passing over the saddle point is studied in the 
field of the metastable potential. The barrier escaping probability is found to be greatly related 
to the fractional exponent a. Properties are revealed to move reversely in the opposite direction 
of diffusion when a is relatively large despite of the zero-approximating effective friction of the 
system. This is very anomalous to the standard Brownian motion. 
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I. INTRODUCTION 


Fractional Brownian motion (fBm) is a quintessential model for the stochastic processes 
which is characterized by long-memory stationary and Gauss-distributed increments Uft 
Different from the standard Brownian motion (sBm), fBm is generally believed to originate 
from the fractional Gaussian noise (fGn) [ 4 , || whose autocorrelation is characterized by 
the Hurst exponent H. For which, H = 1/2 corresponds to the standard Brownian motion, 
H < 1/2 and H > 1/2 denote the sub- or super-diffusive cases, respectively. Theoretically, 
no matter fBm or sBm, both of them can be categorized in the Langevin unification of 
fractional motions [b] in the study of anomalous diffusions |7-f|. 

Since it offers an alternative model of random processes displaying the Joseph effect 
fBm and its variants have attracted considerable attention in the past few decades 


10 ], 


11 ] 


Enormous efforts have been made on the establishment and development of relevant theories. 
However, comprehensively review on the previous studies one can End, little information has 
been got concerning on the dynamical details of a fBm particle diffusing passing over the 
saddle point of a matastable potential energy surface (PES). This maybe because of the 
fractional calculus difficulties immersed in the procedure of solving the kinetic differential 
equations. But recently we hnd that this can be easily achieved from another road of 
mathematical deriving. Therefore in this paper, we report our recent study on this point. 

The paper is organized as follows: in Sec Jill kinematic relations of the diffusion particles 
are obtained by Laplacian solving the fractional Langevin equation. The probability of 
successfully escaping from the potential well is computed in Sec lIIIl for a characteristic 
visualization. Sec II VI serves as a short summary of present results where some discussions 
are also made for a further consideration. 

II. FRACTIONAL LANGEVIN EQUATION AND ITS SOLUTION 


Let us begin from an a 
Langevin equation (GLE) 


ternative approach to fBm basing on the fractional generalized 
12 - 




d 2 


-2H 


= F(x) - ^r(2H - 1 + £*(*), 


dt 2 ~ 2H ' 


( 1 ) 
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where £ H (t) is the fractional Gaussian noise (fGn) which is zero mean with autocorrelation 
(f H (*i)f H (f 2 )> = 2D h H(2H - l)|ti - t 2 \ 2H ~ 2 . D h = [T(l - 2H)cos(H7r)]/(2Hn) is the 
diffusion coefficient identified by the Gamma function T(^) = ( 0 °° t z ~ 1 e~ t dt. 

Although Eq. fTTTl is the standard presentation of fBm, it is mathematically intractable 
due to the fractional calculus. But fortunately we find in our recent study that it can 
be transformed to a more conventional form of fractional GLE by using of the Gamma 
function T(s) and the relation between the external force F{x) and the potential energy 
U(x). Mathematically it can also be derived from a generalized system-plus-reservoir model 


of harmonic oscillators [1 

d 2 x(t) 

and more conveniently 


i.e. 


= F(x) — 7 


dx(t) 


{t-t') a dt' 


dt' + £(£) 


( 2 ) 


mx + / rj(t — t')x{t')dt' + d x ll(x) = £(t) 


(3) 


where r)(t) = 7] a t~ a / r(l — a) is the frictional kernel with fractional exponent 0 < a < 1 
and rj a the strength constant. £(i) is the renewed fGn and (£(£)£(£')) = — t'\~ a the 

fluctuation-dissipation theorem 16]. 

The dominant advantage of Eq. ((2]) lies in its easy processing. The solution of it, namely 
also the equation of motion for the diffusing particle, can be easily obtained by a series of 
Laplacian transformation. For an example, in the particular case of an inverse harmonic 
potential U(x) = —^mco 2 x 2 , we find after some algebra 


x(t ) = (x(t)) + / H(t- r)^(r)dr, (4) 

J o 

in which the mean position of the particle along the transport direction is given by 

(x(t)) = [l + u 2 [ H(T)dr}x 0 + H(t)vo, (5) 

J o 

and the variance of x{t) reads 

cr 2 x (t) = f dtiHit-F) f dt 2 (£(ti)€(t 2 ))H(t-t 2 ), (6) 

Jo Jo 

where H(t) = £ _1 [(s 2 + srj(s) — oo 2 )^ 1 ] is namely the response function and T](s) = r] a s a ~ l 
is the Laplacian transformation of the friction kernel r)(t). 
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III. DIFFUSION DYNAMICS OF FRACTIONAL DAMPING SYSTEM 


In the studies here and following, we extend the results aforesaid to the more realistic 
case of a matestable potential to reveal the diffusion dynamics of a fractional damping (FD) 
system described by Eq. (J 


U(x) = 


17]. This means to set 

Ub — \mujl(x — Xb) 2 near barrier; 

near well. 


(7) 


| mu>Z(x - x a ) 2 

where Ub denotes the barrier height and oji [i = a, b ) indicates the frequency near the 
potential well as well as the saddle point. Although some new derivations should be made, 
but noticing that U(x) can be approximated to be an inverse harmonic potential in the 
neighbourhood of the saddle point, this would not be much of a problem. 

For the study of diffusion dynamics, one of the basic task is to compute the successful 
rate of the particle escaping from the potential well. Mathematically, it always emerges to 
be a complementary error function due to the Gaussian property of noise £(t) and linearity 
of the GLE. i.e. 

(x(t)) 


POO 

P(t)= / W(x,t)d. 

Jo 


x = -erfc 


yielding from a short integration on the reduced distribution function W(x, t) where 


( 8 ) 


W(x,t) = 


a/2tt A u (t) 


exp 


(x 


(x(t))f 


2 A n (f) 


(9) 


is resulted from the joint probability density function (PDF) of the system 18] 


W(x,v,t) = 


3 -|(y/ t (*) A 


2vr|A(t)| 1 / 2 

with y(t) the vector [x — (x(t)),v — (u(t))] and A (t) the matrix of second moments 

- 4 n(i) = al(t) = <[i- {i(«))] 2 >, 

Ai2(t) = An(t) = <[l - (l(t))][K - («(*))]), 

A22O) = al(t) = <[w- (v(t))] 2 ). 


( 10 ) 


(11a) 

(lib) 

(11c) 


Now everything is ready before tracking the trajectories to End the dynamical details. 

In the calculations here and following, we rescale all the quantities so that dimensionless 
unit such as k B T = 1.0 is used. Firstly in Figs. Q] and |2] we plot the time dependent 
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FIG. 1: Time dependent varying of (. x(t )) at various a. Parameters in use are rj a = 2.0, oj = 
ksT = 1.0 , xq = —1.0 and vq = 2.0. 



FIG. 2: Time dependent varying of An(t) at various a. 

varying of (x(t)) and Au[t) (namely also at various a. From which we can see 

that ( x{t)) increases quickly when a is relatively small. But as the increasing of a, (x(t )) 
suddenly decays into a negative diverging function of a. The critical point lies in the nearest 
neighbourhood of a = 0.55 in the particular case of what is considered here. This is a very 
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nontrivial result, because from the view point of PDF evolution, (x(t)) indicates the center 
of the wave packets and dn(t) the width. Therefore it reveals that the center of the wave 
packet may move in the opposite direction of diffusion. 

However, the value of An(t) is always positive increasing. This means, as the forwarding 
of the center the width the wave packet is always expanding. Therefore, despite of the 
unusual movement of the center one may always expect a steady barrier escaping probability. 
As is shown in Fig. [3] the probability for a particle to diffusing passing the saddle point 
tends to be a steady one (defined as P st ) in the long time limit no matter what is the value 
of a. The occurrence of such a nontrivial result may probably be caused by the property of 
long-range correlations (the “Joseph effect”) immersed in the standard process of fBm [2], 



FIG. 3: Time dependent varying of P(t) at various a. 

In order to get more details about the fBm, we plot in Fig. [4] the time dependent varying 
of P st (the steady value of P(t)) as a function of the fractional exponent a. From which 
we can see that P st decays monotonously as the increasing of a. At first thought, this is 
considered to be caused by a strong resistance from the friction kernel. However, when 
we check through the time-dependent varying of r]{t), a startling result is found that the 
effective friction of the system tends to be zero as the increasing of a. As is shown in Fig. 
[5] for more explicitness. But what is the real reason then? Zero damping will lead to low 
barrier escaping, right? 
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FIG. 4: Stationary value P st plotted as a function of a. 



FIG. 5: Friction kernel rj{t) plotted as a function of t at various a. 

In order to answer these questions, we make a thorough investigation on the kinematical 
properties of the system. We find that the mean velocity of the particle can be obtained 
from a short derivation from (x(t)), i.e. 

(v(t)) = ^(x(t)> = H(t)x o + j t H(t)v 0 . (12) 

And succeeding calculation reveals that (v(t)} tends to be positive when a is relatively small, 
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but will turn to be negative at large a. As is shown in Fig. [ 6 ] for more details. The critical 
point is identical to what we have found hereinbefore. Therefore we have known the reasons 
now: for the fBm process with large a, the effective friction of the system is very strong 
at the beginning. Particle inserted in this dissipative environment loses its initial kinetic 
energy quickly. After a very short period of time, it gets a reversed kinetic energy to move in 
the opposite direction of diffusion despite of the zero-approximating damping. This results 
in a low probability of barrier escaping (small P st ) in the long time. 



FIG. 6 : Friction kernel 77 (f) plotted as a function of t at various a. 


IV. SUMMARY AND DISCUSSION 


In conclusion, we have studied in this paper the diffusion dynamical properties of the frac¬ 
tional damping system which may induce fBm. Several anomalous phenomena are revealed 
such as the reverse moving of the particles at large a, the low Pst at zero-approximating 
effective friction and so on. All these results have shown us a vivid picture much different 
from the standard Brownian motion. 

Since theoretically fBm may be a model particularly relevant to subcellular transport, 


unbiased translocation 
solutions 


19 


20 ], the dispersion of apoferritin proteins in crowded dextran 


21] and lipid molecules in lipid bilayer 22]. Therefore we believe that deep research 


will be greatly indeed in the near future in some disciplinary held such as sub-diffusing mRNA 


















molecules 


23], RNA proteins and chromosomal loci within E. coli cells 24], And the results 


we have got here will in no doubt motivate the continuous demystification of these seemingly 
simple subjects. 
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